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ABSTRACT 

Massive black holes in active galaxies are immersed in huge concentrations of late type stars 
in the galactic bulges and also early type massive stars in the nuclear stellar clusters which are 
additionally surrounded by quasi-spherical several kpc scale halos containing from a few hundred 
up to several thousand globular clusters (GCs). It is expected that significant numbers of red 
giant stars, massive stars and also GCs can move through the jet expelled from the central 
engine of active galaxy. We consider collisions of stars from the galactic bulge, nuclear cluster 
and globular clusters with the jet plasma. As a result of such collisions, multiple shocks are 
expected to appear in the jet around these compact objects. Therefore, the plasma in the kpc 
scale jet can be significantly disturbed. We show that particles can be accelerated on these shocks 
up to the multi-TeV energies. TeV leptons emit synchrotron radiation, extending up to the X- 
ray energies, and also comptonize radiation produced in a stellar cluster and also the Microwave 
Background Radiation to TeV 7 -ray energies. We show that such non-thermal radiation is likely 
to be detectable from the intermediate scale jets of the nearby active galaxies for reasonable 
number of stars and GCs immersed within the jet. As an example, we calculate the expected 
non-thermal emission in the X-ray and gamma-ray energies from the nearby radio galaxy Cen A 
from which a steady gamma-ray emission with the complex spectrum has been recently reported 
by the Fermi and the HESS Observatories. 

Subject headings: Galaxies: active — galaxies: individual: Centaurus A — globular clusters: general — 
stars: massive — radiation mechanisms: non-thermal — gamma-rays: general 


1. Introduction 

The non-thermal X-ray emission from large 
scale jets in nearby radio galaxies of the Fanaroff- 
Riley class I is expected to be common feature 
of such type of objects starting from its discov¬ 
ery in two close galaxies Cen A or M 87 (Feigel- 
son et al. 1981, Biretta et al. 1991). This X-ray 
emission smoothly extends up to the radio energy 
range supporting the idea on its common origin in 
the synchrotron process of TeV electrons (Hard- 
castle et al. 2001). TeV electrons are expected to 
meet strong soft radiation field, produced either 
in the synchrotron process or by the stars in the 
galactic bulge and also the Microwave Background 
Radiation (MBR), fro its comptonization to TeV 


energies (e.g. Stawarz et al. 2003, Hardcastle & 
Croston 2011). Observations of persistent TeV 7 - 
ray emission from such sources should provide im¬ 
portant constraints on the high energy particles in 
the kpc scale jets. However, since such electrons 
has to lose energy on a short time scale, they re¬ 
quire an in situ acceleration process. Therefore, 
efficient acceleration of electrons should occur not 
only within the inner jet but also in the kpc scale 
jets. In fact, the observed X-ray knots in the in¬ 
termediate scale jets seems to provide conditions 
for electron acceleration to TeV energies. These 
might correspond to shocks formed in collisions of 
winds from compact objects with the jet plasma 
(e.g. Hardcastle et al. 2003). 

A jet, launched from the vicinity of the cen- 
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tral super-massive black hole (SMBH), has to 
pass through the stellar bulge and globular clus¬ 
ter (GC) halo around central engine. Since the 
number of stars in the central stellar cluster is 
expected to be huge, many of them have to be 
immersed in the jet plasma. In fact, it has been 
argued that as a result of the star-jet collisions, 
stellar winds can provide large amount of barion 
matter into the jet (e.g. Komissarov 1994, Bow¬ 
man et al. 1996, Laing & Bridle 2002, Hubbard 
& Blackman 2006). Shocks formed in collisions of 
stellar winds with the plasma in the inner jet can 
be also responsible for the acceleration of electrons 
which, comptonizing stellar radiation, might be re¬ 
sponsible for the production of 7 -ray flares (e.g. 
Bednarek & Protheroe 1997, Barkov et al. 2010, 
Bosch-Ramon et al. 2012, Araudo et al. 2013, 
Wykes et al. 2014). Also collisions of relativis¬ 
tic jet plasma with dense clouds can turn to the 
7 -ray production (e.g. Blandford & Kdnigl 1979, 
Dar & Laor 1997, Beall & Bednarek 1999, Pur- 
mohammad & Samimi 2001, Araudo et al. 2010, 
Barkov et al. 2012). These clouds might be frag¬ 
ments of the close to the jet supernova explosion 
Blandford & Konigl (1979). Large number of ob¬ 
jects, moving in the galactic bulge, are expected 
to be immersed within the kpc scale jets. 

Galaxies of different types are surrounded by 
spherical halos of GCs with the radii of several 
kpcs. GCs are expected to produce winds which 
are the products of the mixture of the relativis¬ 
tic winds from a large population of millisecond 
pulsars (MSPs) within the GC and slow barionic 
winds from the red giant stellar population within 
GCs (Bednarek & Sobczak 2014). The mixed mil¬ 
lisecond pulsar and stellar winds, expelled from 
the GC, interact with a relatively rare relativistic 
intermediate scale jet. As a result, a double shock 
structure forms around GC. We argue that parti¬ 
cles can be accelerated on such shocks to multi- 
TeV energies. They produce non-thermal radia¬ 
tion by interacting with the magnetic field in the 
jet (the synchrotron radiation) and with the soft 
radiation from the galactic bulge at the sub-kpc 
distance scale and also with the MBR at the a 
kpc distance scale via inverse Compton process. 

In this paper we explore the hypothesis that rel¬ 
ativistic electrons appear in the intermediate scale 
jet as a result of their acceleration on the shocks 
formed in collisions of many compact objects (red 


giants, massive stars or GCs) with the jet. It is 
assumed that at the kpc scale distances jets are 
already non-relativistic and the Doppler boosting 
effects are small. Under this assumption, we cal¬ 
culate the synchrotron and inverse Compton 7 - 
ray emission produced in the intermediate jet for 
different parameters of such scenario. As an ex¬ 
ample, we show that the high energy TeV 7 -ray 
component in the Cen A spectrum, observed by 
the HESS Collaboration (Aharonian et al. 2009, 
Sahakyan et al. 2013), might originate within the 
intermediate scale jet as a result of acceleration 
of electrons in collisions of many compact objects 
with the jet plasma. 

2. Stars and their clusters around AGN 
jets 

Active galaxies are complicated systems sur¬ 
rounded at larger scale by a few important ele¬ 
ments. The central SMBH, with the mass in the 
range ~ 10 ® — 10 ^° Mq, is surrounded by the bulge, 
which has quasi-spherical shape, containing late 
type stars within the radius of the order of ~kpc. 
The mass of the bulge is related to the mass of the 
SMBH and its luminosity is of the order of a few 
tens percent of the whole galaxy. Many SMBHs 
are surrounded by central clusters of luminous and 
massive stars. Such clusters have usually extended 
shape with the main axis in the galactic plane. Its 
mass and dimension is supposed to be determined 
by the destruction process of another galaxy by 
the massive galaxy harboring SMBH. Finally, ac¬ 
tive galaxy is surrounded by a halo of GCs within 
the distance scale of several kpc. The number of 
GCs is related to the mass of the parent galaxy and 
the mass of the SMBH. The parameters of these 
stellar clusters around nucleus of active galaxy are 
discussed below. 

Nuclear clusters of massive stars are observed 
not only around SMBHs in massive active galax¬ 
ies but also in the center of our Galaxy. In the 
case of our Galaxy such cluster is relatively small. 
It has a half mass radius of ~ pc and the mass 
of ~ 2.5 X 10^ Mq (e.g. Schodel et al. 2014) In 
the case of the nearby active radio galaxy Cen A, 
the massive stellar cluster has been formed around 
50-100 Myrs ago as a result of a merger process 
with a gas reach galaxy. The observed star for¬ 
mation region is centered on the Cen A nucleus. 
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It has the visual dimensions of ^ 8x3 kpc^ (Neff 
et al. 2015). The current star formation rate has 
been estimated on 2 M 0 yr“^ which will result in 
production of (6 — 12 ) x 10 ^ Mq of young stars 
(Wykes et al. 2014). Assuming the standard ini¬ 
tial mass function of stars, we estimate the total 
number of stars with masses above 20 Mq in the 
nuclear cluster in Cen A on ~ 3 x 10®. A small part 
of these stars can pass through the jet expelled 
from the super-massive black hole. The rough es¬ 
timate of the number of massive O and WR type 
stars, contained within the 1.5 kpc from the base 
of the jet with the opening angle of the order of 
0.1 rad for the cluster with dimension mentioned 
above, is of the order of a hundred. This simple es¬ 
timate shows that many massive stars can interact 
with the jet plasma. 

Central regions of galaxies are also surrounded 
by the quasi-spherical concentrations of late type 
stars (galactic bulges) with typical dimensions 
of the order of ^1 kpc and luminosity of sev¬ 
eral percent of the total luminosity of the par¬ 
ent galaxy. Interestingly, the masses of bulges 
are related to masses of the central black holes 
(for review on scaling formulas see e.g. Czerny & 
Nikolajuk 2010). In the case of Cen A, the num¬ 
ber of bulge stars within the jet has been estimated 
on ~ 8 X 10® (Wykes et al. 2014). If only a small 
number of these stars are in the red giant phase, 
of the order of ~ 10 “®, then the total number of 
red giants in Cen A jet can be as large as ~ 10®. 

As we mentioned in the Introduction also GCs 
are expected to enter the jet from time to time 
at the kpc distances from the base of the jet. In 
fact the number of observed (or estimated) GCs 
around galaxies seems large enough to guaran¬ 
tee encounters with the jet. For example, in the 
case of our Galaxy more than 150 has been dis¬ 
covered (Harris 2003). Larger galaxies are ex¬ 
pected to have more GCs. The Andromeda galaxy 
is expected to be surrounded by a halo of ~500 
GCs (Barmby & Huchra 2001) but the giant ra¬ 
dio galaxy, M87, contains already ~13000 GCs 
(McLaughlin et al. 1994). It is shown that in ac¬ 
tive galaxies, containing super-massive black holes 
(SMBHs), the number of GCs in the galactic halo 
is correlated with the mass of the central black 
hole (e.g. Spitler & Forbes 2009, Burkert & 
Tremaine 2010, Harris et al. 2014). For example, 
in the radio galaxy Cen A the number of GCs is 


estimated on ^ 1550 (Giiltekin et al. 2009). 

3. Interaction of compact objects with 
AGN jets 

We consider active galaxy which has a jet prop¬ 
agating from the SMBH. The jet has a conical 
structure on the intermediate scale. The jet open¬ 
ing angle at a kpc distance scale is expected to 
be of the order of 0 ~ 0 . 1 d_i rad, see e.g. ob¬ 
servations of the nearby active galaxy Cen A (e.g. 
Clarke et al. 1992, Hardcastle & Croston 2011). 
Then, the solid angle subtended by the jet is, 
AQ = 2.5 X Since the number of com¬ 

pact objects surrounding central engine can be 
very large (see above), many of the late type stars, 
massive stars, and GCs can pass through the jet 
(see Fig.l on the left). As we mentioned above the 
number of these compact objects is shown to be re¬ 
lated to the mass of the central black hole. There¬ 
fore, it is expected that very massive black holes, 
which are on average more luminous, are also sur¬ 
rounded by a larger number of stars and stellar 
clusters. For example, we estimate the number of 
GCs within the jets of two well known nearby ac¬ 
tive galaxies, Cen A and M87, on ~ 10 and ~70 
GCs for the jet opening angle equal to 0.15 rad. 

The power of jets in the nearby radio galaxies 
are typically of the order of ^ 10 ^® erg s“® (e.g. 
in Cen A, Wykes et al. 2013). The ram pres¬ 
sure of the jet plasma can be estimated from (e.g. 
Bednarek & Protheroe 1997), 

^1 = Jy. - 1-1 X (1) 

where the power of the jet is Lj = 10 ^®L 43 erg s“^, 
the distance from the base of the jet is I = II 3 kpc 
and c is the velocity of light. Compact objects, 
mentioned above, are also surrounded by winds 
which pressure can be balanced by the pressure 
of the jet plasma. As a result, a double shock 
structure is expected to appear around the com¬ 
pact objects within the jet as already discussed in 
the case of different types of stars. In the case of 
GCs, hadronic winds from significant population 
of red giant stars mix efficiently with the winds 
from many millisecond pulsars (MSPs) within the 
GC. As a result, the wind emanates from the GG 
with the power determined by the energy output 
from the MSP population within the GG. This 
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Fig. 1.— Schematic representation (not to scale) of the interaction of the winds from compact objects (red 
giants, massive stars, globular clusters) with the plasma of a relativistic jet in active galactic nuclei. On the 
left: large number of stars and GCs from halos around the central engine on the distance scale up to ^kpc 
and several kpcs, respectively. Some of these objects enter the jet expelled by the central super-massive black 
hole. On the right: The example case of collision of a GO with a jet. GCs are expected to contain up to a 
few hundred of millisecond pulsars (MSPs) and a few hundred thousand of normal stars, many of them in 
the red giant phase. The MSP winds and red giant winds mix efficiently within the cluster producing mixed 
pulsar-stellar wind from the GC. This GC wind interacts with the jet plasma. As a result, a double shock 
structure is formed. Leptons can be accelerated at the shock up to multi-TeV energies. These electrons lose 
energy on the synchrotron process and on the IC scattering of the soft radiation produced by stars in the 
bulge (at kpc distance scales) and on the Microwave Background Radiation during their propagation along 
the kpc scale jet. 
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mixture of red giant/MSP winds collide with the 
jet’s plasma as shown schematically in Fig. 1 (on 
the right). We argue that particles can be acceler¬ 
ated at the shocks from the jet site to multi-TeV 
energies. They produce non-thermal radiation by 
interacting with the magnetic field in the jet (the 
synchrotron radiation) and with the soft radiation 
from the galactic bulge at the kpc distance scale 
and also with the Microwave Background Radia¬ 
tion (MBR) at the distance scale up to a few tens 
of kpcs (the inverse Compton radiation). Below 
we consider separately the case of collision of a 
star with the jet and a GC with the jet since these 
types of objects differ significantly in the parame¬ 
ters describing their winds. 

3.1. Stars in the galactic bulge and nuclear 
stellar cluster 

As we argue above two different classes of stars 
are expected to surround the central engine of ac¬ 
tive galaxy at distances up to the order of a kpc. 
They are massive stars from the nuclear stellar 
cluster and late type stars (red giants) from the 
galactic bulge. In fact, collisions of the winds from 
massive stars with the inner jet plasma in active 
galaxies, as a possible scenario for acceleration of 
particles and their radiation, have been already 
discussed in Bednarek & Protheroe (1997). More 
recently, distraction process of red giant stars as a 
result of collisions with the inner jets has been also 
considered as the mechanism for y-ray flares (e.g. 
Barkov et al. 2010 and Bosch-Ramon et al. 2012). 
Such destruction process is possible only in the 
case of very powerful jets and at a few pc scale dis¬ 
tances from the base of the jet. At large distances, 
where the pressure of the jet plasma is weak, the 
red giant star winds are strong enough to stop the 
jet plasma above the red giant surface. As a result, 
the scenario for collisions of red giant stars be¬ 
comes similar to the scenario of collisions of mas¬ 
sive stars with jets on the intermediate scale dis¬ 
tances from the base of the jet. Recently, collision 
of different types of stars with the intermediate 
scale jets, as the mechanism for the non-thermal 
synchrotron emission from jets, has been discussed 
by Wykes et al. (2014) with the application to Cen 
A. 

First, we consider collisions of stars with pa¬ 
rameters of the red giant stars. The scenario for 
collisions of massive stars will be in principle the 


same. We will comment on possible differences at 
the end of this section. We assume that a small 
part of the bulge stars, of the order of ~ 10 “^, 
are in the phase of the red giant. As an exam¬ 
ple, we consider red giants with following typical 
parameters, the radius Rrg = 5 x 10^^ cm, the 
surface temperature Trg = 3700 K, the luminosity 
Lrg = 10^ L©, the mass loss rate Mw = 10 “'^M _7 
M© yr“^, and the wind velocity Uw = 3 x 10®V 3 
cm s“^. As discussed in Sect. 2 , a large number 
of red giant stars, of the order of A^rg = lO^iVg, is 
expected to be immersed in the intermediate jet of 
active galaxy on the distance scale of ~ 1 kpc (see 
estimates for the numbers of bulge stars within the 
jet in Wykes et al. 2014). 

These red giant stars produce stellar winds 
which pressure can be estimated from. 


p* _ 

^ ~ 47ri?2 


1.6 X 10-i®M_7V3Rr^ (2) 


where the distance from the star is i? = li?i pc. 
The stellar wind pressure balances the pressure 
of the jet plasma at the distance from the stellar 
surface at (Bednarek & Protheroe 1997), 


« 3.6 X 10^®(M_7V3)^/^6»_iZ 3/A43^ cm. (3) 


For small distances from the base of the jet, the 
pressure of the jet might not be balanced by the 
pressure of the wind of the red giant above its 
stellar surface. In fact, the shock is above the red 
giant surface for the stars which enter the jet at 
the distances from its base which are larger than, 

I « O.15Lli^/[0.iiM_rV3y^^] pc. (4) 


For the jet powers as observed in nearby radio 
galaxies (e.g. Cen A), L = lO"*® erg s“^, this dis¬ 
tance scale is only of the order of I ~ 0.1 pc. It 
is much smaller than the radius of the galactic 
bulge. Therefore, most of the red giants entering 
the jet in Cen A from its galactic bulge will be 
surrounded by the shocks clearly at some distance 
from stellar surface. The atmospheres of these red 
giants will not suffer direct collisions with the jet 
plasma as considered by Araudo et al. (2010) and 
Barkov et al. (2010). Instead, the jet plasma will 
be stopped at some distance from the red giant 
surface creating a shock as considered by Bednarek 
& Protheroe (1997) for the case of massive stars 
entering the inner jet in active galaxies. 
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The maximum power, that can be extracted 
from the jet by such a shock structure around a 
single star, is estimated on, 

4*h = i^j(^)"«l-45xl035M_7^3 (5) 

Note that this power depends only on the param¬ 
eters of the red giants star. Since the total num¬ 
ber of red giants in the jet can be quite large, a 
significant energy of the jet can be extracted and 
eventually transferred to relativistic particles in 
the acceleration mechanism operating at the colli¬ 
sion region of the jet and stellar winds. 

The parameters of the massive stars, which can 
also enter the jet, are clearly much more extreme 
than expected for the red giants. For example, the 
terminal velocities of the winds of massive O type 
stars can be ?;w ~ (1 — 3) x 10^ km s“^ and the 
mass loss rates Mw ~ 10 “® M 0 yr“^ and of the 
WR type stars ~ (1 — 5) x 10^ km s“^ and the mass 
loss rates ~ ( 0.8 — 8 ) x 10 “^ Mq yr“^, respectively 
(Lang 1992). Stars in different phases of evolution 
can have very large mass loss rates of the order 
of 10“® — 10“^ Mq yr“^, e.g. Asymptotic Giant 
Branch stars or Luminous Blue Variable stars (see 
Table 1 in Wykes et al. 2014). These massive stars, 
although more rare than red giants, will produce 
shocks in jets which are able to extract ~3-4 or¬ 
ders of magnitude more energy from the jet plasma 
than a single red giant. 

3.2. Globular clusters in the galactic halo 

As described in Sect. 2, GCs are expected to 
create a halo on several kpc distance scale around 
the central engine of active galaxy. A substantial 
number should be submerged in the jet with the 
opening angle of the order of ~0.1-0.15 rad. GCs 
are composed from old low mass stars which signif¬ 
icant number, on average about ^ 100 , is expected 
to be in the red giant phase. These red giants 
are expected to be the main contributors of the 
mass into the GC through their stellar winds. In 
fact, the mass loss rate of isolated red giant is ex¬ 
pected to be in the range between 10 “® Mq yr“^ 
to 3 X 10“^ Mq yr“^ (Boyer et al. 2008, Meszaros 
et al. 2009). The mass loss rate in the range of 
10“^ Mq yr“^ to 3 X 10“^ Mq yr“^ is estimated 
from 100 red giants in GC. It is expected that the 
single solar mass star should lose about ~ 0.3 Mq 
during its evolution path to the White Dwarf (e.g. 


see discussion in Heyl et al. 2015). Simple esti¬ 
mates of the mass loss rates in a single GC, with 
the total number of stars equal to 10 ®, gives us 
the average mass loss rate during 10 Gyr equal to 
~ 3 X 10“® Mq yr“^, consistently with the above 
estimated value. In our example calculations, we 
use the mass loss rates of stars in GCs close to 
this value. On the other hand, GCs contain many 
millisecond pulsars which has been accelerated to 
short periods due to the transfer of momentum 
from the stellar companions. These MSPs have 
typical surface magnetic fields of the order of ~ 10 ® 
G and rotational periods of the order of ~3 ms. 
The number of such pulsars in a single GC can 
be of the order of ^ 100 , as estimated from their 
observed 7 -ray luminosity (Abdo et al. 2010c). 
These MSPs lose rotational energy in the form of 
the relativistic pulsar winds that should mix ef¬ 
ficiently with slowly moving winds from the red 
giants (Bednarek & Sobczak 2014). As a result, a 
barion loaded pulsar wind slows down to the sub- 
relativistic velocities. The velocity of the mixed 
pulsar-stellar wind is estimated to be. 


,,GC 


I 2Lmsp 
Mqc 


5.6x 10"(^)i/2 


cm 

I 

s 


( 6 ) 


where the mass loss rate of stars within the GC 
scales with Mac = 3 x 10“®M_5 Mq yr“^ and 
Lmsp = 3 X 10®®L36 erg 8“^ is the rotational en¬ 
ergy loss rate of all millisecond pulsars within the 
GC. 

The pressure of the GC wind is estimated from 
(Bednarek & Sobczak 2014), 

pGC ^ _ 9.5 X 10~^°(L36M_5)^/^ erg 

‘iirB? R\ cm® 

The location of this shock in the jet around the 
GC, i?sh, can be estimated by comparing the pres¬ 
sure of the wind from the GC with the ram pres¬ 
sure of the jet plasma, 

« O.93(L36M_5)i/^0-iZ3/A4f PC- (8) 


The radius of the shock should be larger than the 
core radius of the GC. Note that, most of GCs 
(~ 90%) in our Galaxy have the core radius below 
1 pc (Harris 1996). Therefore, for the parameters 
considered in this paper we expect that this condi¬ 
tion is usually fulfilled for the GCs at the distance 
larger than ~1 kpc from the central engine. 
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The shock structure, formed in collision of the 
jet with GC wind, is expected to have a double 
structure with different conditions on both sites 
(see Fig. 1 on the right). The shock from the jet 
site is semi-relativistic (the velocity of plasma after 
the strong shock drops by a factor of three) and 
has relatively strong magnetic field. On the other 
hand, the shock in the GC wind is non-relativistic 
with weak magnetic field. We expect that only 
shock in the jet plasma can accelerate efficiently 
particles to large energies transferring significant 
part of the jet power to relativistic electrons. The 
power which can be extracted by the shock from 
the jet plasma can be estimated on, 

4°'^ = ^j(^)'~8-6xl038(L36M_5)i/2 (9) 

Note that the power estimated above for the sin¬ 
gle collision of the GC with the jet is clearly larger 
than in the case of a single collision of the red gi¬ 
ant star with the jet but it becomes of the order of 
the power expected in the collision of a single mas¬ 
sive star with the jet. Therefore, collisions of GCs 
with jets at the distance scale of a few kpcs from 
the central engine can be also responsible for the 
formation of powerful shocks in jets at such large 
distances. In the next section we discuss possible 
acceleration of electrons on the shocks formed in 
collisions of compact objects with jet plasma in a 
relatively nearby active galaxies. 

4. Acceleration of electrons at the shock 
in the jet 

In order to estimate the parameters of particles 
which could be eventually accelerated in such sce¬ 
nario, we consider the physical conditions within 
the jet. The upper limit on the magnetic field at 
the base of the jet can be estimated by assuming 
that the Poynting flux (Tp) through the jet, is a 
part (/i) of the total jet power (Tj), i.e. 

Lp = TrrfjjC/BcT^ = fiLj = erg s“^, (10) 


field at the base of the jet can be estimated by 
reversing Eq. m, 

Bb ~ 520(MT43)^/V(rM8) Gs. (11) 


The perpendicular component of the magnetic 
field in the jet (in respect to the jet axis) at the 
distance, I, from its base can be approximated by, 


B (0 


Bhnn 

(nil + ^0 


1.7 X lO"'^ 


T9.ih 


Gs, (12) 


for large distances from the base of the jet. In 
this case, B{1) does not depend on the mass of the 
black hole. 

We assume that electrons are accelerated at the 
shock from the site of the jet plasma to energies 
limited by their synchrotron energy losses. The 
energy gain rate by electrons is assumed to be de¬ 
termined by the so called acceleration coefficient, 

e = 

Bacc = ^cE/Rp « 9.1 X 10®C-3S(0 eV s-\ (13) 


where E = IBxeV TeV is the energy of electrons 
and i?L is the Larmor radius of electrons. This ac¬ 
celeration process can be saturated by the energy 
losses of electrons. The synchrotron energy loss 
rate is. 


Esyn — 7 :CO'tUb{1){ - 

3 TTleC^ 


^1/43 Big V eV 
(r0_ii3)^ si ’ 


where the energy density of magnetic field is, 


Ub{1) = BV(87r) R. 7A0^iL43/{Te_,kf — , (15) 


ctt is the Thomson cross section and me is the 
electron rest mass. Then, the maximum energy 
of electrons, limited by the synchrotron energy 
losses, are 


K, 


55(4)^/^ Rs 135 
B 


(e-3^-iZ3r)i/^ 


TeV. (16) 


where the inner radius of the jet can be related 
to the SMBH mass Hn = 3rsch = lOi^^Mg cm, 
the mass of the BH is Mbh = lO^Mg, the energy 
density of the magnetic field is Ub = B^/Stt, and 
r is the Lorentz factor of the jet. The jet can be 
either Poynting flux dominated (then /i ~ 1) or 
matter dominated (then /i ^ 1). The magnetic 


Electrons with these maximum energies cannot 
have the Larmor radii larger than the dimension of 
the shock around the compact objects. The con¬ 
ditions, i?L = and Rp — introduce an¬ 
other absolute limit on the maximum electron en¬ 
ergy, equal to E^^ « 1.4 x /T 

TeV and El k 1.8 x /T TeV. The 
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limits due to the Larmor radius are clearly above 
the maximum energies of accelerated electrons. 
However, another limit on the maximum energies 
of electrons is introduced by the time scale spend 
by electrons at the shock region. This limit differ 
for stars and GC. Therefore, we consider it sepa¬ 
rately. 


4.1. Stars in jet 

The maximum energies of electrons can be also 
limited by their advection time scale along the 
shock with the flowing plasma. This time scale 
can be estimated from, 


"^adv — 


3R 


sh 


/3c 


s,(17) 


where /3 is the velocity of the jet plasma in units of 
the velocity of light. The comparison of the accel¬ 
eration time scale (race = E/Eg,cc, see Eq.[T3|) with 
the advection time scale, Tadvj gives us the limit 
on the maximum energies of electrons accelerated 
in the stellar shock, 

EaV ~ 5.6C_3(/iM_7n3)'/V(/?r) TeV. (18) 


Note that advection limit on the energy of elec¬ 
trons mainly depends on the parameters of the 
stellar wind. On the other hand, the synchrotron 
limit mainly depends on the parameters of the jet. 
In general, the advection limit becomes important 
in respect to the synchrotron limit for the red gi¬ 
ant stars in contrast to the case of the massive 
stars. 

Electrons accelerated at the shock around the 
star can also lose efficiently energy on the in¬ 
teraction with the stellar radiation. In order 
to evaluate importance of this process we cal¬ 
culate the energy density of stellar photons at 
the shock region, [/* = asBT^(R*/Rth)'^ ~ 

1.7 X 10 ^TiRl^Li3/[M-7V3ie.ih)^] eV cm’^, 
where the surface temperature of the red giant 
type star is T* = 3.7 x 10 ^T 3 K and its radius 
i?* = 5 X 10 ^^i?i 2 cm, and ctsb is the Stefan- 
Boltzmann constant. On the other hand, energy 
density of the magnetic field at the shock from the 
jet site is given by Eq. (15). The relative energy 
losses of electrons on these two processes (in the 
Thomson regime, i.e. for Eg -C 28 O/T 3 GeV) de¬ 
pends only on the ratio of the energy densities of 
the radiation and the magnetic field. Therefore, 


electrons with energies clearly below ^ 300 GeV 
will lose energy more efficiently on the IC process 
than on the synchrotron process. However, the 
energy losses of electrons becomes comparable for 
TeV energies to which these particles are expected 
to be accelerated at the shocks provided that the 
jet is Poynting flux dominated (/r ~ 1). Since the 
winds of massive stars are clearly more extreme 
than those of red giants, a relatively stronger IG 
7 -ray emission (in respect to the synchrotron emis¬ 
sion) is expected in the case of massive stars within 
the jets. 

The power in relativistic electrons extracted 
from the shock, created in collision of the jet with 
stellar winds, is estimated on, 

L: = « 1.45 x (19) 

where t] = 0 . 1 ry_i is the energy conversion effi¬ 
ciency from the shock to relativistic electrons, and 
iV* = lO^iVe is the number of stars within the jet. 
The energy conversion efficiency, 7 , is usually ex¬ 
pected to be in the range ^0.1-0.2. The number of 
massive stars within the jet is expected to be a few 
orders of magnitudes lower than the number of red 
giants (see Sect. 2). However the product of the 
mass loss rate and the wind velocity of a massive 
star is typically 3-4 orders of magnitudes larger 
(see Sect. 3.1). As a result, the power in relativis¬ 
tic electrons accelerated on the shocks formed by 
the whole population of red giants and massive 
stars is expected to be on the similar level. 

4.2. Globular clusters in jet 

As above, we constrain the maximum energies 
of electrons due to their advection along the shock 
around the GG within the jet. We compare the ac¬ 
celeration time scale, race, with Tadv = 

This condition limits energies of electrons to, 

E^g ^ 430^_3M'/'(i36M_5)'/V(/3r) TeV. (20) 

In general, this limit is less restrictive than the 
advection limit for stars within the jet, since 
the shock dimensions around GGs are expected 
to be clearly larger for typical parameters of 
the scenario. The synchrotron limit on the ac¬ 
celeration of electrons dominates over the ad¬ 
vection limit for the following condition on the 
least constrained parameter of the model, i.e. 
the acceleration coefficient should be, ^_3 > 




0.16»_iZ3/3?rV(i43i36M_5)i/VM^/^- If the above 
condition is fulfilled, then we can estimate the 
characteristic energies of synchrotron photons pro¬ 
duced by electrons on Ex ~ WeC^(i3/i?cr)(£^syn/WeC^)^ ~ 
1.4 X 10®^ keV, where i?cr = 4.4 x 10^^ Gs is the 
critical magnetic field strength. In the case of 
Cen A, the X-ray emission up to ^ 10 keV is 
observed (Hardcastle et al. 2006). Therefore, the 
acceleration coefficient is expected to be at least 
^ > 10“^. Electrons with such energies should 
produce synchrotron emission extending to X-ray 
energies but also TeV y-rays by scattering optical 
radiation produced within the galactic bulge and 
the MBR. 

Let us check whether electrons, accelerated at 
the shock around GC, will manage to cool effi¬ 
ciently on the synchrotron or the Inverse Comp¬ 
ton (IC) processes. The electron cooling mean free 
path on the synchrotron process in the magnetic 
field of the jet is, 
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( 21 ) 


The synchrotron cooling process of electrons will 
occur locally in the jet, provided that the syn¬ 
chrotron mean free path becomes shorter than 
the characteristic distance scale which we iden¬ 
tify with the distance of the injection place of 
electrons from the base of the jet, i.e. Agyn < I- 
The above condition is fulfilled for electrons with 
energies E > {Li^T) TeV. The IC en¬ 

ergy losses of electrons can in principle be deter¬ 
mined by the soft radiation produced by the GC 
itself, the galactic bulge thermal radiation and 
the MBR. The bulge radiation might play an im¬ 
portant role in the inner part of the intermediate 
scale jet (i.e ~kpc scale). The MBR can start to 
be important in the outer parts of the large scale 
jet when the energy density of the magnetic field 
drops significantly. The energy density of radi¬ 
ation from the GC at the location of the shock 
can be estimated on C/qc = I^Gc/47rc(i?^^)^ ~ 
32OL6L43/(L36M_5)^/^(0_il3)^ eV cm“^, where 
the optical luminosity of the GC is Lgc = 
= 4 X 10^® erg s“^ and Lq is the 
luminosity of the Sun. The energy density of 
the bulge radiation is estimated on C/buige = 
ibuige/(47ri?g^^jggC) Ri 70Lii/i?§ eV cm”^^ where 
the bulge luminosity is Lbuige = Lq, and 

the bulge radius is i?buige = I-R 3 kpc- The energy 


density of the MBR is C/mbr = 0.25 eV cm“^. 

The mean free path for electron energy losses on 
IC process in the Thomson regime in these radia¬ 
tion fields can be estimated on, 

^Ic = cE /Eic ~ 94/(ETevGev) kpc, (22) 

where C/rad = IGev eV cm“®. We expect that 
multi-TeV electrons will mainly scatter stellar ra¬ 
diation in the Klein-Nishina regime but the esti¬ 
mate based on the Thomson cross section allows 
us to have impression about importance of the IC 
energy losses. The IC losses becomes important 
when the mean free path of electrons is compara¬ 
ble to the characteristic distance scale on which 
they propagate. In the case of the radiation field 
from the GC, is always clearly lower than XJq, 
for reasonable parameters of the scenario. There¬ 
fore, energy losses of electrons in the GC radiation 
can be safely neglected. In the case of the bulge 
radiation, the mean free path of the IC energy 
losses can become comparable to the character¬ 
istic size of the bulge (of the order of kpc) and 
also to the mean free path of electrons on the syn¬ 
chrotron process. Therefore, the IC energy losses 
in the bulge radiation should be taken into account 
on the kpc distance scale from the base of the jet. 

The mean free path for the IC energy losses of 
electrons in the MBR is of the order of ~kpcs, for 
electrons with energies of a hundreds TeV. There¬ 
fore, IC losses of electrons are expected to become 
important in the outer parts of the jets at the dis¬ 
tance scale of several kpc from the base of the jet. 

Note that the GC halo still extends to such dis¬ 
tances. 

We conclude that in fact, electrons can be ac¬ 
celerated to large energies at the shocks formed in 
collisions of GCs with jet plasma. They should 
produce synchrotron radiation extending up to X- 
ray energy range and also multi-TeV 7 -rays by 
scattering radiation from the galactic bulge and 
the MBR. 

The total power in relativistic electrons, accel¬ 
erated at the shocks around GCs, can be estimated 
on, 

LfC = rjNGcLg^ - 8.6 x (23) 

where the shock power from a single GC is given 
by Eq. O and Vqc = lOTVi is the number of GCs 
within the jet. 
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The maximum energies of electrons are limited 
in this case by their synchrotron energy losses to 
values given by Eq. (16). For the expected values 
of the parameters describing considered scenario, 
the power in electrons, accelerated at the shocks 
around GCs, is estimated to be about an order 
of magnitude lower than in the case of the shocks 
around stars. However, electrons, accelerated at 
the shocks around GCs, are injected into the jet 
at much larger distances from its base (i.e. above 
~kpc). At these distances the jet magnetic field 
strength is already weak. Therefore, electrons lose 
relatively more efficiently energy on the IC process 
(producing 7 -rays) than on the synchrotron pro¬ 
cess (producing radiation below X-ray energies). 

In the next section we perform detailed calcu¬ 
lations of the multi-wavelength spectra expected 
from the intermediate scale jet assuming that dif¬ 
ferent types of compact objects (red giants, mas¬ 
sive stars, GCs) form multiple shocks in the jet 
plasma. As an example, the parameters of the 
nearby radio galaxy, Gen A, are applied. 

5. Non-thermal radiation from electrons 

We calculate expected synchrotron and IC spec¬ 
tra produced by electrons accelerated on the fronts 
of multiple shocks within the jet for different pa¬ 
rameters describing compact objects and the jet 
content. It is assumed that electrons are acceler¬ 
ated with the differential power law spectra above 
some minimum energy Amin. The spectral index 
equal to —2 is selected for the example calcula¬ 
tions as expected in the shock acceleration sce¬ 
nario. The maximum energies of electrons, Amax, 
are determined either by the synchrotron energy 
losses (Eq. (THl) or by the escape from the accelera¬ 
tion region due to the advection of electrons with 
the jet plasma (Eq.[T 8 ]or[^. These electron spec¬ 
tra have been normalized to the power transferred 
to electrons from the shocks in the jet (see Eqs.[T9] 
and E51) . In order to obtain the synchrotron and 
the IG 7 -ray spectra, we inject electrons at differ¬ 
ent distances from the base of the jet and simulate 
their propagation in the jet magnetic field and the 
radiation fields, discussed above, by applying the 
Monte Garlo method. The Klein-Nishina effect 
has been taken into account when calculating the 
IG 7 -ray spectra. In these example calculations, 
we assume that the jet is semi-relativistic, with the 


apparent velocity equal to /3app = 0.6 and moves 
at relatively large angle to the observer’s line of 
sight estimated on a = 50° in the case of Gen A. 
For the above parameters the velocity of the jet is 
estimated on /3 0.52 and its Lorentz factor on 

r « 1.172. The Doppler factor of the jet in Gen 
A at the distance of a hundred pc from its base is 
then estimated on A = l/[r(l —/3cosa)] Rs 1.282. 
Note that at the kpc distance scales the jet may 
significantly decelerate. Therefore, Doppler effects 
can be negligable in the example case of Gen A. 
Such assumption clearly simplifies the calculations 
of the non-thermal radiation presented in this pa¬ 
per. The jet has a constant opening angle equal 
to 9 = O.I rad and the jet power is Aj = 10"^^ erg 
s“^. For these basic parameters, we investigate the 
spectra as a function of two other parameters de¬ 
termining the relativistic electrons in the jet, the 
jet magnetization parameter /i and the accelera¬ 
tion efficiency of electrons in the jet These two 
parameters determine the relative part of energy 
of relativistic electrons lost on the synchrotron and 
on the IC process. 

5.1. Emission from star-jet collisions 

The calculations of the IC 7 -ray spectra in the 
case of collisions of stars with the inner jet (i.e. at 
distances < 1 pc) have been already considered by 
e.g. Bednarek & Protheroe (1997). In such case, 
the shock is located relatively close to the massive 
companion star and relativistic electrons initiate 
cascade in the dense stellar radiation field. The 
calculations of the synchrotron emission from elec¬ 
trons accelerated on the multiple shock structures 
around the stars in the intermediate scale jets has 
been recently considered by Wykes et al. (2014). 
But no IG 7 -ray emission has been calculated in 
this model. 

We consider the star-jet interaction process on 
a large distance scale from the base of the jet. In 
such case, the cascade process initiated by elec¬ 
trons accelerated around massive stars is ineffi¬ 
cient since the shock is too far from the star. 
Electrons escape with significant energies from the 
shock and propagate mainly in the jet volume. At 
this stage, electrons lose energy on the synchrotron 
process in the magnetic field of the jet and on the 
IG process by scattering different radiation fields 
such as radiation of compact star, diluted radia¬ 
tion from the galactic bulge and the MBR. Elec- 
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Fig. 2.— The synchrotron and IC spectra (SED - spectral energy distribution) produced by electrons 
accelerated at the multiple shocks in the jet plasma around the red giants (on the left) and massive stars (on 
the right), which entered the jet homogeneously on the range of distances from the base of the jet, Lmin = 10 
pc and Lmax = 1 kpc. The results are shown for the Poynting dominated jet, for the value of magnetization 
parameter equal to /x = 1 (dotted curves)) and the matter dominated jets described by ^ = 0.1 (dashed), 
and /X = 0.01 (solid). The upper panel shows the spectra for the acceleration parameter of electrons equal 
to ^ = 10“^ and the bottom panel for ^ = 10“^. The parameters of the jet are the following Lj = 10"^^ erg 
sft 6* = 0.1 rad. It is assumed that 10® red giants and 100 massive stars are present in the jet at the same 
time. The energy conversion efficiency from the jet to relativistic leptons is 10%. The other parameters of 
the stars and the radiation field are reported in the main text. 
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trons propagating outside the bulge, with typical 
dimension of i?buige = 1 kpc, will see the radiation 
field diluted by the factor (V^buige)^- In these 
example calculations we assume that electrons are 
accelerated on shocks around stars entering uni¬ 
formly the jet between 10 pc up to 1 kpc. Note 
that in such case maximum energies of electrons 
injected from the shocks into the jet depend on the 
distance from the base of the jet as given by Ea.fTHl 
or Eq. [THl The parameters of stars have been fixed 
on Mw = 10“^ Mq yr“^ and u* = 30 km s“^ (for 
the red giant stars) and on Mw = 10 “^ M© yr“^ 
and Vi, = 10 ^ km s“^ (for the massive stars). 

Results of these example calculations of the 
synchrotron and IC 7 -ray spectra are shown in 
Eig. 2. We investigate the spectra for the Poynting 
flux dominated jets (the magnetization parameter 
of the jet /r = 1 ) and for the matter dominated jets 
(described by ^ = 0.1 and 0.01). We also show the 
spectra obtained for different values of the accel¬ 
eration efficiency of electrons ^ = 10 “^ and 10 “^. 
In all considered cases, the synchrotron spectra 
dominate over the IC spectra. As expected, the 
IC 7 -ray spectra are stronger for less magnetized 
jets. In the case of red giants, the cut-off in syn¬ 
chrotron spectrum is determined by the advection 
time scale of electrons along the shock structure. 
In contrast, the synchrotron spectra cuts-off in¬ 
dependently on the jet magnetization for shocks 
around massive stars. 

5.2. Emission from globular cluster-jet 
collisions 

Collisions of GCs with jets in active galaxies 
have not been considered up to now as a source 
of disturbances in the jet plasma (shocks) which 
might be responsible for acceleration of particles. 
However, a large number of GCs around active 
galaxies, and evidences that their numbers are cor¬ 
related with the mass of the central black hole 
(Harris et al. 2014), suggest that significant num¬ 
ber of GCs is immersed in the jet. We showed 
above that shocks around GCs provide conditions 
for acceleration of electrons to hundreds of TeV. 
Applying the numerical code for the propagation 
of electrons in the jet and their radiation described 
above, we calculate the synchrotron and IC spec¬ 
tra. Since such calculations have not been con¬ 
sidered before, we also show the results of calcu¬ 
lations of the spectra produced at different dis¬ 


tances from the base of the jet. The spectra are 
investigated as a function of parameters describ¬ 
ing the acceleration process of electrons (i.e. the 
magnetization parameter of the jet /i and the ac¬ 
celeration efficiency of electrons ^). Electrons are 
injected at fixed distance from the base of the jet 
and, after leaving the shock, they are frozen in the 
jet plasma. In such case, the synchrotron energy 
losses of electron are determined by the profile of 
the magnetic field along the jet (see Eq. [T^. On 
the other hand, these electrons also comptonize ra¬ 
diation from the galactic bulge (with photon den¬ 
sity dropping with the distance along the jet) and 
the MBR. Note that electrons are accelerated to 
larger energies farther along the jet. As in the case 
of collisions of stars with the jet, we assume that 
electrons reach the power law spectrum extending 
up to the maximum energies described above. In 
the case of GCs, the minimum energy of electrons 
is fixed on 1 TeV which corresponds to the char¬ 
acteristic energy of leptons injected by the mil¬ 
lisecond pulsars into the mixed pulsar/stellar wind 
from the GC. 

In the example calculations of the synchrotron 
and the Inverse Gompton spectra, we apply the pa¬ 
rameters derived for the nearby radio galaxy Gen 
A: the mass of the black hole Mg = 0.5, the jet 
power L 43 = 1 , the jet opening angle 0 -i = 1 , the 
jet velocity /3 = 0.5 (T Ri 1.15 ~ 1), and the pa¬ 
rameters of typical for GGs: GG stellar luminosity 
Lg = 1, the power supplied by MSPs Tag = 3, the 
mass loss rate of red giant stars M _5 = 3, and 
the parameters of the host galaxy bulge, the bulge 
stellar luminosity Tguige = the bulge ra¬ 

dius i ?3 = 1. The acceleration coefficient is taken 
to be ^_3 = 1 in order to be consistent with the ob¬ 
servations of the non-thermal X-ray emission from 
the jet. The radius of the shock around GG is then 
Rsh ~ 2 . 7 Z 3 pc, the magnetic field strength along 
the jet is B{r) Ki 5.6 x /l^, Gs, the maxi¬ 

mum energies of electrons due to synchrotron en¬ 
ergy losses are E^yn ~ 170(Z3//r)^/^ TeV, the max¬ 
imum energies of electrons due to advection from 
the shock are Eadv ~ 1800/r^/^ TeV. The compar¬ 
ison of these two last limits gives us the range of 
distances for which synchrotron limit is more re¬ 
strictive, i.e. I 3 < 100/r^. We conclude that the 
maximum electron energies are always limited by 
synchrotron losses in the case of Poynting domi¬ 
nated jet. They are limited by the advection pro- 
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Fig. 3.— The synchrotron and IC spectra (SED) produced by leptons accelerated at the shock in the jet 
plasma around the GC which enters the jet at different distances from its base, I = 1 kpc (on the left), 
3 kpc (middle), and 10 kpc (on the right). The results are shown for the Pointing dominated jet, for the 
value of magnetization parameter equal to ^ = 1 (dotted curves), and the matter dominated jets for /i = 0.1 
(dashed), and fj, = 0.01 (solid). The upper panel is for the acceleration parameter of electrons equal to 
^ = 10“^ and the bottom panel for ^ = 10“^. The parameters of the jet are as in Fig. 2. It is assumed that 
20 GCs are present in the jet and the energy conversion efficiency from the jet to relativistic leptons is 10%. 
The parameters describing the GC and the radiation field of the bulge are reported in the main text. 
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cess in the case of the matter dominated jet with 
jjL = 0.01. For ^ = 0.1 the acceleration process 
of electrons is limited by the synchrotron energy 
losses in the inner part of the jet, i.e. below I 3 = 1 , 
but in the outer parts of the jet by the advec- 
tion process. The mean free paths of electrons are 
Asyn ~ 240?§/(/xiiiTev) pc On the synchrotron pro¬ 
cess, w 1.3/iilTeV kpc on the IC scattering 

of bulge radiation , and A^®^ « 380/iilTeV kpc 
on the IC scattering of the MBR. The characteris¬ 
tic distance scales for IC energy losses of electrons 
on the border between the Thomson (T) and the 
Klein-Nishina (KN) regimes in the bulge radiation 
are comparable to the bulge dimension. These dis¬ 
tance scales for IC energy losses of electrons with 
the TeV energies in the MBR are comparable to 
the dimensions of jets in AGNs. 

In Fig. 3, we show the synchrotron and the In¬ 
verse Compton spectra for different distances from 
the base of the jet, I = I kpc, 3 kpc, and 10 kpc. 
The Poynting flux dominated jets (/i = 1) and the 
matter dominated jets = 0.1 and 0 . 01 ) are con¬ 
sidered. In the case of the Poynting dominated 
jets, the IC spectra are on a much lower level 
than the synchrotron spectra especially at lower 
distances from the base of the jet. On the other 
hand, the IC spectra start to dominate over syn¬ 
chrotron spectra at larger distances from the base 
of the jet in the case of matter dominated jets de¬ 
fined by the magnetization parameter /r ~ 0 . 01 . 
Therefore, most of the energy of relativistic elec¬ 
trons can be transferred to the TeV 7 -rays in jets 
at distances of a few kpcs from its base. The 
synchrotron spectra, produced in the jet by these 
electrons, clearly extend up to the X-ray energy 
range, even for the advection dominated accelera¬ 
tion process of electrons, for the considered range 
of parameters, 0.01 < ^ < 1 and 10 “^ < C < 10 “^. 

The synchrotron and IC spectra, produced by 
electrons accelerated at shocks formed by GCs en¬ 
tering the jet at the range of distances from its 
base between I kpc and 20 kpc, are shown in 
Fig. 4. It is assumed that GCs enter the jet homo¬ 
geneously over such range of distances and ejects 
electrons with the spectra and powers as described 
above. The results are shown for the parameters 
as in previous figures. It is clear that in the case 
of the matter dominated jets the level of the TeV 
7 -ray emission is comparable to the level of the 
synchrotron X-ray emission. In this case, the 7 - 


ray spectra peak at multi-TeV energies and the 
synchrotron spectra extends above ~I0 keV. How¬ 
ever, in the case of Poynting dominated jets, the 
synchrotron emission clearly dominates over the 
TeV 7 -ray emission. The 7 -ray luminosity is ex¬ 
pected to be about an order of magnitude below 
the X-ray luminosity and the 7 -ray spectra peak 
at GeV energies. 

We show that in the case of matter dominated 
jet models, the collisions of GCs with the kpc 
scale jets can be responsible for the detectable syn¬ 
chrotron X-ray emission and also for the TeV 7 - 
ray emission. The relative levels of the emission 
in these two energy ranges should provide impor¬ 
tant constrains on the content of jets at kpc scale 
distances. Note that in the matter dominated jets 
the maximum energies of accelerated electrons are 
determined by their escape from the shock region 
due to the advection process. These electrons are 
not completely cooled when propagating in the jet 
on the distance scale of 20 kpc considered in the 
paper. Therefore, the IC spectra are so hard be¬ 
low the peak at TeV energies. The lower energy 
electrons are cooled on much longer distance scale. 
They should produce steeper spectra on similar 
angular scales in the case of distant AGNs for 
which the contribution to the IC spectrum comes 
from more extended part of the jet. Therefore, 
such hard spectra are only expected from the in¬ 
termediate scale jets in the nearby AGNs such as 
Cen A for which the distance of 20 kpc corresponds 
to angular size of the source several arc min. We 
considered only semi-relativistic jets as observed 
from nearby AGNs. In the case of relativistic jets, 
expected in the distant AGNs, these maximum en¬ 
ergies of accelerated electrons are expected to be 
lower (see Eq. 20). Therefore, the spectra may 
not extend to ~I0 TeV. These spectra are not 
expected to be hard enough to provide interest¬ 
ing constraints on the Extra-galactic Background 
Light. 

6. Application to Centaurus A 

Centaurus A is the closest radio galaxy (dis¬ 
tance 3.8 ±0.1 Mpc) associated with the elliptical 
galaxy NGC 5128 (Harris et al. 2010). The SMBH 
in Cen A has the mass estimated on Mbh = 
(5.5±3.0) X lO”^ M 0 (Cappellari et al. 2009). Two 
jets, clearly asymmetric on kpc scale, are observed 
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Fig. 4.— As in Fig. 3 but for the case of GCs which 
entered the jet at the range of distances, ^min = 1 
kpc and ?niax = 20 kpc. The results are shown 
for the Poynting dominated jet, for the value of 
magnetization parameter equal to /x = 1 (dotted 
curves) and the matter dominated jets for /x = 
0.1 (dashed), and ^ = 0.01 (solid). The spectra 
are calculated for the acceleration parameter of 
electrons equal to ^ = 10 “^ (upper panel) and for 
^ = 10“^ (bottom panel). The other parameters 
of the model are these same as in Fig. 3. 


in Cen A. The viewing angle of the jet has been 
estimated on ^50° (Tingay et al. 1998, Hardcas- 
tle et al. 2003). The projected speed of the jet at 
hundred pc distance scale has been measured on 
0.5c (Hardcastle et al. 2003). The power of the jet 
has been estimated on Pj ~ 10"^^ erg s“^ (Wykes 
et al. 2013). 

Due to its proximity, X-ray emission is clearly 
seen from the kpc scale jet in Cen A (Feigelson et 
al. 1981). At the distance between ~l-3 kpcs, the 
Chandra observations (0.4-2.5 keV energy range) 
show that the diffusive X-ray emission dominates 
over X-ray emission from the knots. This X-ray 
emission has been interpreted as due to the syn¬ 
chrotron process (Hardcastle et al. 2006). Such 
synchrotron radiation requires the presence of the 
TeV electrons in the kpc scale jet. They have to be 
accelerated close to the radiation site due to short 
energy loss time scales. Recently, a point like 7 - 
ray source with the power low spectrum (spectral 
index ~2.7) extending to ^2-3 GeV has been de¬ 
tected by Fermi satellite from the core of Cen A 
(Abdo et al. 2010a). The analysis of the 4 years of 
the Fermi data shows that 7 -ray emission extends 
up to ~ 50 GeV. However, the spectrum flattens 
above ~4 GeV (Sahakyan et al. 2013). The spec¬ 
tral index changes from 2.74 ± 0.03 below 4 GeV 
to 2.09 ± 0.20 at higher energies. This flattening 
of the 7 -ray spectrum is consistent with the detec¬ 
tion of the TeV source towards Cen A by the HESS 
Collaboration (Aharonian et al. 2009). It suggests 
that the TeV and hard GeV emission might origi¬ 
nate in this same mechanism. The lower energy 
points of these TeV observations (> 250 GeV) 
links nicely to the hard high energy component de¬ 
tected by Fermi but the spectral index of the TeV 
emission is again better described by a steep spec¬ 
trum (spectral index 2.73 ± 0.45stat ± 0.2syst, see 
Aharonian et al. 2009). Such curious spectral be¬ 
havior, lack of detected clear variability, and pure 
angular resolution of 7 -ray observations suggest 
involvement of at least two radiation processes or 
emission regions (Sahakyan et al. 2013). Due to 
the large angular extend of the Cen A jets on the 
sky (of the order of ~10°), the Fermi Observa¬ 
tory was also able to discover GeV 7 -ray emission 
from the regions of the giant radio lobes formed in 
collisions of jets with intergalactic medium (Abdo 
et al. 2010b). All these high energy observations 
show that the high energy processes are char- 
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acteristic not only for the direct vicinity of the 
SMBH but also for the large scale jets. We won¬ 
der whether the highest energy component (at ~1 
TeV), observed from Cen A by the HESS Collab¬ 
oration, can originate in the kpc scale jet dne to 
the acceleration of electrons on the shocks in the 
jet which are produced in the collisions of compact 
objects (stars, GCs) with the jet plasma. 

In order to test this hypothesis, we performed 
calculations of the synchrotron and the IC spec¬ 
tra expected in terms of the above jet/compact 
objects collision model, applying the known pa¬ 
rameters of the jet and the surrounding medium 
of Cen A. The results of calculations are compared 
with the observed high energy spectrum from Cen 
A (Fig. 5). In this figure the X-ray spectrum rep¬ 
resents emission from the kpc scale jet (Hardcastle 
et al. 2006). The Fermi spectrum comes from the 
central part of the galaxy harboring Cen A (Sa- 
hakyan et al. 2013). The origin and location of this 
steady spectral component is at present unknown 
due to the pure angular resolution of the Fermi- 
LAT telescope. Finally, the data points show the 
measurements of the TeV y-ray emission from Cen 
A inner/intermediate scale jet. We assumed that 
the red giant stars and massive stars enter the jet 
uniformly at distances between 10 pc and 1 kpc 
from the base of the jet. Electrons are accelerated 
at the shocks with the energy conversion efficiency 
of ?7 = 10%. The spectra are calculated for the 
presence of 10® red giants (dotted curves) and 200 
massive stars (dashed curves) within the jet. The 
example parameters of these stars are described 
in Sect. 5.1. For the considered parameters of 
the acceleration model (^ = 0.01, ^ = 10“®), the 
emission from electrons accelerated on the shock 
around red giants is not expected to contribute 
significantly to the high energy spectrum observed 
from Cen A. On the other hand, about a hundred 
of massive stars in the jet can produce synchrotron 
emission on the level observed from the jet of Cen 
A. The accompanying y-ray emission, from the 
componization of background radiation by these 
same electrons, is on the level of the TeV y-ray 
emission reported by HESS from Cen A. We also 
show the synchrotron and IC spectra produced in 
collisions of 20 GCs with the jet at larger distance 
scale (dot-dashed curves) for these same parame¬ 
ters of the jet. The synchrotron emission is not ex¬ 
pected to contribute significantly to the observed 


X-ray spectrum. However, the IC y-ray emission 
peaks at the TeV energies. The cumulative y-ray 
spectrum from collisions of massive stars and GCs 
with the jet is clearly consistent with the level of 
the TeV y-ray emission observed by HESS. There¬ 
fore, such model might be responsible for the high¬ 
est energy y-ray component observed from Cen A. 
If the jet in Cen A still moves with the velocity of 
the order of ~ 0.5c over kpc distances, then the 
calculated spectra (in Fig. 5) should be enhanced 
by a factor of ~ 2.7. As a result, our estimates 
of the number of compact objects within the jet 
should be reduced by this factor. 

We conclude that the steady TeV y-ray emis¬ 
sion, observed from the central part of the radio 
galaxy Cen A, can originate in the kpc scale jet. 
This emission (or at least a part of this emission) 
should be steady independently on the activity 
state of the inner jet. Since the expected emis¬ 
sion is above the level of sensitivity of the future 
Cherenkov Telescope Array (CTA), the question 
on the existence of the low level persistent emission 
from the intermediate scale jet in Cen A should be 
dehnitively answered in the coming years. 

7. Conclusions 

We propose that the compact object/jet colli¬ 
sion model can be responsible for the non-thermal 
emission from the intermediate scale jets observed 
in close radio galaxies. We consider collisions of 
different type of objects, starting from the red 
giant stars from the galactic bulge, through the 
massive stars from the nuclear central cluster and 
finishing on the globular clusters in the galactic 
halo. Collisions produce multiple shocks in the jet 
on which leptons can be accelerated to TeV en¬ 
ergies. These leptons are responsible for the dif¬ 
fusive synchrotron and inverse Compton radiation 
from the kpc scale jets. Contrary to the inner jet 
collision models, leptons accelerated on the shocks 
do not lose energy close to the shock region (near 
the stars) but they are injected into the jet. They 
are transported with the jet plasma to farther dis¬ 
tances from the base of the jet losing energy on the 
mentioned above radiation processes. Therefore, 
this high energy emission is expected to be steady 
during the time scales of the jet dimension (i.e. 
thousands of years), in contrast to the emission 
expected from collisions of stars with the inner jets 
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Fig. 5.— The comparison of the high energy ob¬ 
servations of Cen A jet with the calculations of 
the synchrotron and IC spectra expected from the 
collisions of GCs and stars with the jet in Cen A. 
The X-ray emission from the kpc scale jet in Cen A 
is taken from Hardcastle et al. (2006), the Fermi 
7 -ray spectral measurements are from Sahakyan 
et al. (2013) and the TeV y-ray spectrum mea¬ 
sured by HESS from Aharonian et al. (2009). The 
emission expected from the interaction of 10® red 
giants with the jet (dotted curves) and 200 mas¬ 
sive stars (dashed curves) is calculated assuming 
that these stars enter the jet on the distance scale 
between 10 pc and 1 kpc. The emission from 20 
CCs entering the jet on the distance scale between 
1 kpc and 20 kpc (dot-dashed curves). The energy 
conversion efficiency from the shock to leptons is 
assumed to be 10%. The sensitivity of the CTA 
in marked by the thin dotted curve. 


(below ~1 pc, e.g. Bednarek & Protheroe 1997). 
Moreover, due to the saturation of the accelera¬ 
tion process of leptons by the synchrotron or the 
advection time scales, the larger energy leptons 
are expected to be accelerated at farther along 
the jet. Therefore, non-thermal emission from jets 
should show clear dependence on the energy with 
the lower energy y-rays produced closer to the base 
of the jet. We predict that the multi-TeV y-rays 
should have tendency to be produced farther from 
the base of the jet. In contrast, such feature may 
not be observed in the case of the synchrotron X- 
ray emission if the acceleration of leptons is sat¬ 
urated by synchrotron energy losses. Only when 
the maximum electron energies are constrained by 
their escape (advection) from the acceleration re¬ 
gion, the maximum energies of synchrotron radi¬ 
ation should show energy dependence on the dis¬ 
tance from the base of the jet. 

We compared the example spectra, obtained in 
terms of such a model, with the observations of 
the intermediate scale jet in Cen A (Fig. 5). It 
is found that in general the observed level of the 
X-ray emission can be explained as a result of 
the synchrotron radiation produced by electrons 
accelerated on the shocks around massive stars. 
On the other hand, the TeV y-ray emission is ex¬ 
pected to be produced by leptons accelerated on 
the shocks around massive stars (below ~1 TeV) 
and by electrons accelerated on the shocks around 
CCs (above ~1 TeV). This hypothesis could be 
tested in the near future with the operation of 
the next generation y-ray telescopes such as CTA. 
CTA is expected to reach ~1 arc min resolution at 
10 TeV and the integral sensitivity of ~ 2 x 
TeV cm“^ s“^ at a few TeV (Acharya et al. 2013). 
With such angular resolution the TeV y-ray emis¬ 
sion from the kpc scale jet in Cen A could be re¬ 
solved and the prediction on the energy dependent 
emission with the distance from the base of the jet 
could be confirmed or disproved. 

We applied the collision model for the radio 
galaxies whose intermediate jets are expected to 
be rather slow and observed at a relatively large 
angle to the line of sight. If the jets still moves 
relativistically at the kpc distance scale, then 
the non-thermal radiation can be significantly 
Doppler boosted (see e.g. recent paper by Bosch- 
Ramon 2015). In such case the collisions of differ¬ 
ent types of compact objects (massive stars, red 
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giants, globular clusters) with the relativistic jets 
could provide natural explanation for the power¬ 
ful non-thermal emission from blazars observed at 
small angle to then line of sight. 
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